Exploring classically chaotic potentials with a matter wave quantum probe 
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We study an experimental setup in which a quantum probe, provided by a quasi-monomode guided 
atom laser, interacts with a static localized attractive potential whose characteristic parameters are 
tunable. In this system, classical mechanics predicts a transition from regular to chaotic behavior as 
a result of the coupling between the different degrees of freedom. Our experimental results display 
a clear signature of this transition. On the basis of extensive numerical simulations, we discuss the 
quantum versus classical physics predictions in this context. This system opens new possibilities for 
investigating quantum scattering, provides a new testing ground for classical and quantum chaos 
and enables to revisit the quantum-classical correspondence. 

PACS numbers: 03.75.Pp,05.45.Mt,37.25.+k,41.85.Ew 



The realization and the manipulation of Bose-Einstein 
condensates (BECs) have led to many experiments on 
atom-atom interactions in the low energy regime: low 
partial waves elastic scattering has been analyzed to de- 
duce accurate determination of collisional properties pi 
2] , the possibility of tailoring the scattering length using 
either magnetic [3[ or optical Feshbach resonances [4J has 
been demonstrated, and the role of the dimensionality on 
atom-atom scattering processes through very anisotropic 
confinement has been studied Q . 

In all these examples, the BEC is used to reveal the 
microscopic atom-atom interactions properties. BECs 
have also been placed into disordered @, Ej or non- 
commensurable bi-chromatic potentials [9] . Such experi- 
ments do not study a single matter wave scattering event. 
This is to be contrasted with experiments where a BEC 
bounces on a mirror [l0[. However, in this case, the po- 
tential with which the BEC interacts is extremely simple. 

In this article, we demonstrate the use of a guided 
matter- wave [lll-KUj to study, in a confined environment, 
the interaction of a matter wave with a static tunable lo- 
calized attractive potential (LAP) associated with a com- 
plex dynamics. The transverse confinement of the guide 
can be well approximated by an harmonic confinement. 
The LAP superimposed to the guide breaks the harmonic 
character of the potential experienced by the atoms and 
couples the longitudinal and transverse degrees of free- 
dom (see Fig. [I]). We study this scattering problem the- 
oretically from the weak to the strong coupling regimes 
for which the system exhibits classical chaos. This can be 
done by tuning the LAP power (strength) and position 
d. Finally, we show that these different regimes can be 
clearly identified in our experimental results. 

Our experiment provides a direct access to the mod- 
ulus squared of the wave function through absorption 
imaging. This is to be contrasted with experiments car- 
ried out with microwave in chaotic structure in which 
the measurement is performed at a well-defined position 




FIG. 1: (color online). Sketch of the experimental setup. 
Atoms outcoupled from a BEC at point A are guided by the 
horizontal beam towards point B. In the vicinity of B, atoms 
also experience the attractive force exerted by the LAP (see 
text) whose center is at a distance d from the guide axis. 



[15j . In mesoscopic physics, experiments exploit the con- 
ductance measurements to infer properties of the chaotic 
structure [l6l-[l9j. However, in these systems, the chaotic 
dynamics effects cannot be singled out from disorder and 
interactions effects. In contrast with pioneering exper- 
iments using cold atoms to investigate classical chaotic 
behavior in static optical billiards [20|, [2l| , in our experi- 
ment we prepared atoms in a nearly pure quantum state 
and we analyzed a single event rather than relying on av- 
erage observables. Other studies on quantum chaos have 
been performed with cold atoms in the different context 
of ID time-dependent lattice potentials [22l424j. 

The experimental setup that delivers guided matter 
waves has been described in detail elsewhere [l3|. A 
rubidium-87 BEC is produced in a crossed dipole trap 
using two focused laser beams with a wavelength of 1070 
nm and waists wq ~40 /im. Atoms are prepared in the 
\F = l,mf = 0) internal state through a spin distilla- 
tion process implemented during the evaporative cooling 
stage [12]. The quantum probe is realized by magnetic 
or all-optical outcoupling from the trap (located at A in 
Fig. [I]) [13]. Atom-atom interactions turn out to be neg- 
ligible in the guided atom lasers produced in this man- 
ner [13[ . Experimentally, the propagating matter wave is 
characterized by three quantities: the atom flux (a few 
10 5 atoms/s), the mean velocity v (tunable from 10 to 
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FIG. 2: (color online), (a) Typical trajectory in the weak cou- 
pling regime (see text). The interaction with the LAP yields 
a downstream propagation with a transverse oscillation, (b) 
Example of trajectory in the strong coupling regime. The 
complex dynamics may lead to the reflexion of the particle. 
In both (a) and (b) the LAP potential is represented by the 
blue shaded area, (c) Example of exponential separation £ of 
nearby trajectories in phase space revealing the chaotic nature 
of the potential. The red dotted curve is the best exponential 
fit. Average over 100 couples of trajectories close to the cen- 
tral one is taken, (d) Amplitude of the transverse excitation 
as a function of the ratio Pl ap / Pg and the displacement d of 
the LAP: transmitted (red colors) and reflected (blue colors) 
particles. A chaotic zone with a complex dynamics can be 
clearly identified. 

30 mm.s -1 ) and the population of the transverse modes. 
Experimental parameters are chosen to have the mean 
excitation number (n) below one, which corresponds to 
a ground state population larger than 50% [l3|. Such a 
probe therefore explores the phase space with a trans- 
verse resolution close to h 2 . 

The horizontal beam used for the dipole trap provides 
the guiding potential which can be approximated by a 
harmonic transverse confinement with a measured fre- 
quency uj±/2tt ranging from 170 Hz to 200 Hz (guide 
power, Pq — 100-150 mW). The LAP originates from 
the same laser source as for the crossed dipole trap and 
is slightly detuned to wash out interferences. It has a 
waist of u>lap ~ 60 /im and an adjustable power, Plap 5 
from to 2W. It intersects the guide beam at a distance 
AB = 700 fim from the trap (A) along the propagation 
direction (see Fig. [1]). 

The interaction of the quantum probe with a well- 
centered LAP (d = 0) has been investigated theoreti- 
cally in [25j. To analyze the experimental results of the 
scattering problem for d ^ it is instructive, as a first 
guideline, to use classical mechanics. 



Consider that atoms are classical point-like particles 
delivered at the trap location (A) and at the bottom of 
the guide (see Fig. [I]). As a result of the longitudinal cur- 
vature of the guide beam, atoms are accelerated towards 
point B. Close to P, atoms feel the combination of the 
guide and the attractive potential. When the LAP is well 
centered {d = 0), it does not generate a transverse oscil- 
lation downstream. However, when the transverse sym- 
metry is broken (d ^ 0), the LAP has a dramatic effect. 
In the weak coupling regime, the beam goes downstream 
and exhibits a transverse dipole oscillation (Fig. (2ji). A 
concomitant reduction of the longitudinal velocity occurs 
due to energy conservation. This output channel is sta- 
ble against small changes of d/w^AP and Plap/^g- In 
the strong coupling regime, the output channel can be ei- 
ther upstream (Fig. [2)3) or downstream depending on the 
different parameters of the LAP (strength and displace- 
ment) and turns out to be extremely sensitive to small 
parameters changes (Fig. [2]i). 

The complex trajectories that result from the strong 
coupling between longitudinal and transverse degrees of 
freedom illustrate the chaotic nature of the potential in 
this regime (see Fig. [2fc) where the Lyapunov exponents 
can be shown to be positive and the distribution of trap- 
ping times is found to exhibit an exponential decay (data 
not shown) . Figure [2]i summarizes the variation of the 
asymptotic transverse amplitude as a function of the ra- 
tio Plap / Pg and the displacement d of the LAP. It ex- 
hibits a region with strong variations of output channels 
with parameters, that we identify as the chaotic zone. 
To further characterize this zone, we have verified that 
the set of boundaries between left and right output chan- 
nels as a function of d displays a fractal structure with 
a fractal dimension £. For a wide variety of fixed ra- 
tio Plap/Pg in the chaotic zone, we find an exponent 
£ « 0.9 revealing a strong fractality of the basin bound- 
aries, usually associated with chaotic scattering [26|. 

However, the incident guided atom laser is not a classi- 
cal probe since it is in the ground state of the transverse 
confinement [131 ] . In addition, the potential experienced 
by the atoms is not harmonic and the dynamics cannot 
thus be exactly mapped onto the classical one. 

To explore the quantum-classical correspondence, we 
have developed a full 2D numerical simulation of the 
wave packet that represents the incident elongated mat- 
ter wave. The extremely large anisotropy of the system 
with a longitudinal length of a few mm and a transverse 
size initially of a few hundreds of nm and the long in- 
teraction time inherent to the complex dynamics makes 
this kind of calculation memory size and time consuming, 
and even unrealistic in 3D with the currently available 
supercomputers. We have thus compared the dynamics 
of 2D quantum wavepackets with simulations of classical 
probability densities of same initial position and momen- 
tum distributions. Parameters are chosen in close simi- 
larity with the experimental conditions. In the weak cou- 
pling regime, the scattering is regular, the quantum wave 
packet keeps its structure and performs coherent oscilla- 




FIG. 3: (Color online), (a) (b) (c) Color density plots of the 
wave function from 2D numerical simulations. Color scale is 
not linear. Point B is the same as in Fig.[T] Wave packets are 
prepared in the ground state of the initial trap, and released 
in the guide at time t = (wg — 45 /im and uj± = 2ty x 203 
Hz) at 700 /im from the LAP (wlap = 41 /mi). The in- 
cident mean velocity of the wave packet is v ~ 10 mm.s -1 
close to B. (a) d = 2 /mi (d/w LA p 0.05), Plap 0.5 W 
(Plap/Pg — 2.78), propagation time t = 120 ms, (b): d = 15 
/im (d/^LAP = 0.37) and Plap = 1W (Plap/Pg = 5.56), 
t 120 ms, (c) d = 20 /im (d/w LA p 0.5), Plap = 1 W 
(Plap/Pg = 5.56), t — 160 ms, (d) projected distribution (c) 
on the z axis of the wave function f dx\ip(x,z)\ 2 (black) and 
the initially identical classical distribution (red/gray); clas- 
sical and quantum curves are identical for cases (a) and (b) 
(data not shown). 

tions downward the LAP (Fig.[3K). In this regime, classi- 
cal and quantum predictions are indistinguishable. In the 
strong coupling regime for which the classical counterpart 
is chaotic, the wave packet loses its filamentary structure. 
In this regime, the numerically simulated wavepacket fol- 
lows the central trajectory for short time (see Fig. [3b). 
For larger time, it loses its structure and classical and 
quantum dynamics can become different (see Fig. [3fc and 
Fig. [3]i). The quantum and classical results differ sig- 
nificantly when the wave function overlaps with itself in 
the LAP region and therefore exhibits some interference 
patterns (see Fig. [3]i). 

The scattering problem depends on many parameters: 
those that characterize the incident wave (velocity, num- 
ber of occupied transverse modes, longitudinal coherence 
length, ...) and those that concern the confinement (pow- 
ers, waists and distance d). Experimentally, we have in- 
vestigated this large parameter space by varying the ratio 
Plap/Pg which in turn slightly modifies the displace- 
ment d (see below). The experiment duration, r ~ 300 
ms, is chosen so that typically half the guided atom laser 
has interacted with the LAP in the weak coupling regime. 
This duration was chosen in order to observe the asymp- 
totic regime in the weak coupling regime and the complex 
dynamics associated with the strong coupling regime. 
Absorption images are taken after time-of-flight once the 
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FIG. 4: (color online). Interaction between a guided atom 
laser and a LAP (experimental results). Absorption image 
(a) is taken in the weak coupling regime and, (b) and (c) 
in the strong coupling regime after a 15 ms time-of-flight. 
(d) Measurement of the displacement d in the weak coupling 
regime: the ratio A/S of the transverse size before and after 
the interaction is plotted vs Plap/Pg (v c± 25 mm.s -1 , and 
uj±/2tt c± 200 Hz). The dot-dashed curve corresponds to 3D 
(classical) simulation with d = 4 /im (d/wLAP ~ 0.07). Most 
experimental points (disks) are inside the pink area delimited 
by the simulated curves (gray lines) for 3.7 and 4.3 /im. 



confinement potentials are switched off. In this way, the 
ratio between the final, A, and initial, 5, transverse sizes 
can readily be extracted. Examples of experimental re- 
sults are shown in Fig. HJ The two regimes (regular and 
chaotic behaviors) are clearly observed when the power 
is increased (see Fig.SH and, Fig.H)3 and c respectively). 
The transverse oscillation observed after the LAP in the 
weak interacting regime corresponds to a coherent trans- 
verse dipole oscillation (Fig. H^i) in close similarity with 
the corresponding quantum simulation (Fig. [3^). The 
apparent decrease of the transverse size of the helix mo- 
tion is due to the guide anisotropy and the projection of 
the trajectory along the imaging axis. 

For given parameters, the displacement d of the LAP is 
reproducible from shot to shot. However, d varies slightly 
with the LAP power. This latter effect is due to thermal 
effect generated by the RF power in the acousto-optic 
modulator that control the LAP power. In the parame- 
ter space of Fig. [2]i, the experimental results thus do not 
follow a vertical line. This is clearly illustrated in Fig|4]i 
where we have represented A/ <$, in the weak coupling 
regime, for both the experimental points, and 3D (classi- 
cal) simulations that assume a fixed displacement d = 4 
jam. Here, we use the fact that classical and quantum 
predictions coincide in this regime. The observed satu- 
ration for Plap > IOPg suggests that in this parameter 
regime d is approximately constant d ~ 4 /im ±0.3 /im. 
For Plap < IOPg, the initial curvature is opposite to 
the simulated one revealing a change of d as Plap is in- 
creased. Using a linear increase of d up to Plap = IOPg 
towards the value d = 4 /im and a constant value after- 
wards gives a very good agreement with the experimental 
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data (see Fig. [Hi). 

Interestingly enough, the experiments turn out to be 
robust with respect to the velocity dispersion of the in- 
coming matter wave. This effect, confirmed by numerical 
simulations, originates from the fact that the larger avail- 
able energy for larger incident velocity is compensated for 
by the smaller interaction time with the LAP. 

Experimentally, we observe two types of images in the 
chaotic regime: those for which the trajectory can be 
well-identified (Fig. HJd) and those with a more diffuse 
cloud (Fig. Ht). Numerical simulations (see Figj3b) sug- 
gest that the former can be attributed to situations for 
which the chaotic dynamics did not have enough time 
to fully develop. The latter may result from either the 
ergodic mixture as in Fig. [3b that occurs in the chaotic 
regime for sufficiently long time or situations in which 
the relatively high local atomic density favors elastic col- 
lisions and thus the possibility of exploring a larger phase 
space zone with various output channels. 

In Fig. [5j we assess the reproducibility of the experi- 
ment by using different runs with same initial conditions, 
where small initial differences are inevitably present. Ex- 
amples of the different experimental images are shown in 
the top inset. We define the quantity <Jij(z) as the dif- 
ference of the mean transverse size of the wave packet at 
coordinate z downward the LAP (located at zb) for two 
different runs i and j taken in the same conditions. We 
note (<j(z)) = N^T 1 ^2 iK j &i,j(z) where N r is the number 
of runs (A^ r = 5 in our experiment) and ((c)) is (cr(z)) 
averaged over z > zb- In Fig. [5] (top), this latter quan- 
tity is represented vs Plap / Pg • When the system enters 
the chaotic (strong coupling) regime a large dispersion 
in the transverse size is observed as a result of the high 
sensitivity of the system to parameter changes and ini- 
tial conditions. The transition to chaos becomes visible 
for Plap/Pg °° 8, compatible with values of d « 7.5 /im 
(d/wLAP = 0.125) (see Fig. EH). 

To make more visible the chaotic regime in the exper- 
iments, in Fig. [5] (bottom) we show the average distance 
(a(z)). As long as the wavepacket moves along the z 
coordinate, this coordinate gives the time evolution of 
the quantum wave function. Unavoidably small differ- 
ences were present in the initial conditions. In the reg- 
ular regime, these differences remain approximately of 
the same order of magnitude, while in the chaotic regime 
they display a sharp exponential amplification when the 
wavepacket enters the chaotic potential. This is rem- 
iniscent of what happens for classical trajectories (see 
Fig. [2fc). Strictly speaking, (cr(z)) does not give a direct 
access to the Lyapunov exponent since more runs would 
be required and we have only a projection of the density 
probability in position space (instead of the phase space 
trajectory in classical physics). However, according to 
our numerical studies, (cr(z)) in many cases can be used 
to give a lower bound to the Lyapunov exponent. 

In conclusion, we have explored the effects of a LAP 
on a guided atom laser. A classical analysis enables to 
identify different regimes, from regular to chaotic and un- 
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FIG. 5: (color online) Experimental difference between quan- 
tum evolutions for different runs with same initial conditions 
(same parameters as in Fig. (Top) ((cr)) (see text) vs 
Plap/Pg- The onset of the chaotic regime is clearly visi- 
ble on both the large value of ((cr)) and its error bar. In- 
set: Examples of these different runs in the weak (left) and 
strong (right) coupling regimes. (Bottom) (<j(z)) (see text) 
vs z. Solid line is Flap = 0.06 W (Plap/Pg = 0.6), dot- 
ted line is Plap = 0.6 W (Plap/Pg = 6), dashed line is 
Plap = 1 W (Plap/Pg = 10). We observe the amplification 
of small initial deviations when the beam enters the chaotic 
zone (z rsj zb ~ 0) as in Fig. [2t. 



stable dynamics. These different regimes are also visible 
in the quantum simulations as well as in the experiments. 
Differences between classical and quantum dynamics are 
especially noticeable in the non-regular regime. This sys- 
tem enables to investigate both classical and quantum 
chaos in a well-controlled environment. Two situations 
should be especially adapted to magnified quantum fea- 
tures: (i) the 2D geometry which favors interferences of 
the wavefunction with itself and (ii) situations for which 
LAP size is on the order of the transverse size of the ini- 
tial wave function. In this latter regime quantum diffrac- 
tion and/or tunneling effects should play an important 
role. The experimental exploration of these two regimes 
provides a natural continuation of our work. Another 
direction for future work includes the study of the non- 
linear influence of atom-atom interaction and in particu- 
lar matter wave solitons [27|, [28| . We think our approach 
paves the way for the use of guided matter waves to study 
highly non-equilibrium dynamics and to use such systems 
as quantum simulators for complex quantum scattering 
problems. 

We thank R. Mathevet and T. Lahaye for useful com- 
ments. We acknowledge financial support from the 
Agence Nationale de la Recherche (GALOP project), 
the Region Midi-Pyrenees, the university Paul Sabatier 



5 



(OMASYC project) and the Institut Universitaire de computers. 
France. We thank CalMip for the use of their super- 



[1] Ch. Buggle, et al Phys. Rev. Lett. 93, 173202 (2004). 

[2] N. R. Thomas, et al Phys. Rev. Lett. 93, 173201 (2004). 

[3] S. Inouye, et al, Nature 392, 151 (1998). 

[4] C. P. Koch, F. Masnou-Seeuws, and R. Kosloff, Phys. 
Rev. Lett. 94, 193001 (2005); R. Ciurylo, E. Tiesinga, 
and P. S. Julienne, Phys. Rev. A 71, 030701(R) 
(2005); G. Thalhammer, et al, ibid 71, 033403 (2005); 
K. Enomoto, et al Phys. Rev. Lett. 101, 203201 (2008). 

[5] D. S. Petrov and G. V. Shlyapnikov, Phys. Rev. A, 64, 
012706 (2001). 

[6] T. Bergeman, M. G. Moore, and M. Olshanii, Phys. Rev. 

Lett., 91, 163201 (2003). 
[7] J. Billy, et al, Nature 453, 891 (2008). 
[8] D. Dries, et al Phys. Rev. A 82, 033603 (2010). 
[9] G. Roati, et al, Nature 453, 895 (2008). 
[10] K. Bongs, et al, Phys. Rev. Lett. 83, 3577 (1999); 
A. S. Arnold, C. MacCormick, and M. G. Boshier, Phys. 
Rev. A 65, 031601 (2002). 
[11] W. Guerin, et al, Phys. Rev. Lett. 97, 200402 (2006). 
[12] A. Couvert, et al, Europhys. Lett. 83, 50001 (2008). 
[13] G. L. Gattobigio, et al, Phys. Rev. A 80, 041605(R) 
(2009). 



[14] G. Kleine Burning, et al, Appl. Phys. B 100, 117 (2010). 

[15] W. Lu, et al Phys. Rev. E 61, 3652 (2000). 

[16] F. Sols, M. Macucci, U. Ravaioli and K. Hess, Appl. Phys. 
Lett. 54, 350 (1989). 

[17] R. A. Jalabert, H. U. Baranger, and A. D. Stone, Phys. 
Rev. Lett. 65, 2442 (1990). 

[18] C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

[19] M. A. Topinka et al, Nature 410, 183 (2001). 

[20] N. Friedman, et al, Phys. Rev. Lett. 86, 1518 (2001). 

[21] V. Milner, et al, Phys. Rev. Lett. 86, 1514 (2001). 

[22] F. L. Moore, et al, Phys. Rev. Lett. 75, 4598 (1995). 

[23] P. Szriftgiser, et al, Phys. Rev. Lett. 89, 224101 (2002). 

[24] C. Ryu, et al, Phys. Rev. Lett. 96, 160403 (2006). 

[25] G. L. Gattobigio, et al New J. Phys. 12, 085013 (2010). 

[26] E. Ott, Chaos in dynamical systems (Cambridge Univer- 
sity Press, Cambridge, 1993). 

[27] L. Khaykovich,e£ al Science 296, 1290 (2002); K. E. 
Strecker, et al Nature 417, 150 (2002). 

[28] J. Martin, B. Georgeot and D. L. Shepelyansky, Phys. 
Rev. Lett. 101, 074102 (2008). 



